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The energy gap is calculated for the ground state quantum computer circuit, which was recently 
proposed by Mizel et.al. When implementing a quantum algorithm by Hamiltonians containing 
only pairwise interaction, the inverse of energy gap 1/A is proportional to A/"**", where A'^ is the 
number of bits involved in the problem, and A'^'' is the number of control operations performed in a 
standard quantum paradigm. Besides suppressing decoherence due to the energy gap, in polynomial 
time ground state quantum computer can finish the quantum algorithms that are supposed to be 
implemented by standard quantum computer in polynomial time. 
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Quantum computer is widely believed to outperform 
its classical counterpart for some classically difficult 
problems I3. Althou gh m an y schemes have been 
suggestedll llTali, H Elll the main obstacle 
is decoherence that makes very difficult the realization 
of even one qubit or one gate for quantum computing. 
Ground state quantum computer (GSQC) is a new ap- 
proach proposed by Mizel et.al.0, [3 which mim- 
ics time evolution of a system by space distribution of 
ground state wavefunction. The advantage of GSQC is 
that appreciable energy gap suppresses decoherence from 
environment as long as temperature is low enough. 

To analyze the performance of GSQC, the key is to 
evaluate the scale of A, the energy gap between the 
ground state |V'o) and the first excited state which 
is related to the time cost. Mizel et.al. gave estimations 
on single qubit and two-qubit circuit in |l4j . However, 
they didn't extend to arbitrary circuit with boost or pro- 
jection Hamiltonians applied, while boost and projection 
Hamiltonians are necessary to make sure that lipo) con- 
centrates on the position corresponding to the final time 
in standard paradigm so that measurement on GSQC can 
read out desired information with appreciable probabil- 
ity. In ^5 it was proposed to shorten length, for a single 
qubit, by inserting teleportation circuits, however, it will 
be shown that there the energy gap scale is incorrect. In 
the present paper, the scaling of energy gap is analyzed 
for general circuit that can be used for any known quan- 
tum algorithm. I will use boundary condition different 
from [13, and show that there exist dangerous cases, in 
which the energy gap shrinks exponentially as the num- 
ber of qubits increases. In order to prevent such expo- 
nential shrink, teleportation circuits are also inserted for 
multiple interacting qubits. Finally energy gap scaling 
for general quantum algorithm is presented. 

A standard computer is characterized by time depen- 
dent state as: \tp{ti)) — Ui\ip{ti-i)), where ti denotes 
instance of the z-th step, and Ui represents for unitary 
transformation. For GSQC, the time sequence is mim- 
icked by the space distribution of the ground state wave- 
function \tpo)- 

As proposed by Mizel et.al. [l^, a single qubit may be 
a column of quantum dots with multiple rows, and each 



row contains a pair of quantum dots. State |0) or |1) is 
represented by finding electron in one of the two dots. 
It is important to notice that only one electron exists 
in a qubit. A GSQC is made up by circuit of multi- 
ple interacting qubits, whose ground state is determined 
by the summation of single qubit unitary transforma- 
tion Hamiltonian h{Uj), two-qubit interacting Hamilto- 
nian h{CNOT), boost Hamiltonian h{B,X) and projec- 
tion Hamiltonian h(\'y),X). 

The single qubit unitary transformation Hamiltonian 
has the form: 



- (cjUjCj^i + h.c. 



(1) 



where e defines the energy scale of all Hamiltonians, 

Cj = c] c] ]^ , cj g is the electron creation operator 

on row j at position 0, and Uj is two dimension matrix 
representing for unitary transformation from row j — 1 
to row j. The boost Hamiltonian is: 



h\B,\) 
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h.c. 
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which amplifies the wavefunction amplitude by large 
number A compared with previous row in jV'o)- The pro- 
jection Hamiltonian is 



h={\l),\) = 
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where I7) represents for state to be projected to on row j 
and to be amplified by A. The interaction between qubit 
a and /3 can be represented by h{CNOT): 



a,j,l^a,j,l'l'/3 
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where for ^ ^, its subscription a represents for qubit a, 
b for the number of row, c for the state |c). With only 
h^{Uj) and h{j^{CNOT), its ground state is;i4]: 

= [l + (l + C| 

M (l + C'l^NCp^,^,) 



(5) 



IS are positr 

nite, and are the same as those in [lia ll4L . Only pair- 
wise interaction is considered for Hamiltonians involving 
multiple qubits. 

The input states are determined by the boundary con- 
ditions applied upon the first rows of all qubits, which can 
be Hamiltonian h'^ = E{I + aiai) with ai being Pauli 
matrix and J^i = 1- For example, if /i° = E{I + (Jz), 
then iV'o) on the first row is if h'^ = E{I — cr^;), then 
it is (|0) -I- |1)) /%/2. Unlike in (Ij, we set £: > e > 0, 
thus boundary Hamiltonians are not perturbation. If E 
is large enough, the energy gap is independent of its mag- 
nitude. A typical value may be -E « lOe. 

Although one can analytically obtain energy gap for 
single qubit with n rows_14j, it's difficult to calculate for 
a complicated circuit. However, we can still manage to 
find its scale. 

Let's first consider the simplest case: a n-row single 
qubit without projection or boost Hamiltonian. The 
Hamiltonian X]"=i ^{Ui) is in fact the kinetic energy just 
like a particle in a one-dimension box with length n, thus 
A oc e/n^. The unique ground state j-^o) is uniform over 
all rows, and \il>i) is orthogonal to \iI>q) on each row. Due 
to large value of i? in ft,*', 1-01) on the first row is nearly 
zero, and it is an mono-increasing function that reaches 
maximum on the final row so as to keep wavefunction 
smooth. When a boost Hamiltonian h{B,\) is applied 
to the qubit's final row, concentrates there, where 
the wavefunction amplitude is about A times those at 
other rows on average, so the wavefunction amplitude on 
other rows is 0(l/-\/A^ +n — \). When the qubit length 
n <C A, iV'i) is nearly a linear function of position ex- 
cept for at final row, and its amplitude at the first row is 
zero. The local kinetic energy is almost constant along 
the whole qubit, proportional to the square of difference 
of wavefunction between neighboring rows, hence energy 
gap scales as: 



A2- 



(6) 



If instead of boost Hamiltonian, a projection Hamilto- 
nian /i(|7). A) is applied on the last row of a single qubit, 
then IV'o) on the final row is restricted to state I7). As- 
suming besides amplitude and phase |V'o) on the second 
last row is at state |^), which is normalized, and (^|7) 
is appreciable, then will concentrate on the final 
row, and should have little weight there, otherwise 



fa) 




FIG. 1: Part (a) shows two qubits interacting with each other 
by h{CNOT). Part (b) shows a chain of qubits interacting 
with each other through h{CNOT). The label / stands for 
identity transformation, CNOT for h{CNOT), U for arbi- 
trary unitary transformation and B for boost Hamiltonian 
h{B,\). 



(V-'i|'0o) 7^ 0. Thus A remains at about e/n^, independent 
of A. 

Numerical calculations on single qubit in various situa- 
tions have confirmed all the above analysis. For example, 
a 6-row qubit with h{B, X) on final row and with bound- 
ary Hamiltonian, = 10e(/ — tiz), on its first row, has 
energy gap A = 0.0782, 0.0174, 1.94 x 10^, 1.96 x 
10-4, 1.96 X lO-^e at A = 1, %/lO, 10, %/lOOO, 100 re- 
spectively. 

Next we consider two qubits interacting with each 
other through h{CNOT) as shown in Fig.l^), with the 
left qubit as control qubit, the right one as target qubit, 
and both ending with h{B,X). By observing IiJjq) in 
Eq.lO, it's easy to see that it has the form 



j control \ 
^upstream/ 



IV', 



control \\ \ I target \ 

downstream/ J \ y^upstreami 



_\ I control x\ . target ^ 

I downstream 1 1 t downstream/ 
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where the Hamiltonian h{CNOT) divides both qubits 
into upstream and downstream parts with h(B, A) in 
downstream parts. h(B, A) on the final row of the target 
qubit also raises the amplitude of \ipo) on the downstream 
rows of the control qubit, while h(B, A) on control qubit 
doesn't influence the amplitude on the target qubit at all 
because both parts of states entangled with the down- 
stream state of control qubit. Thus on the control qubit 
the amplitude of upstream over its final row is l/A^, and 
it is 1/ A on target qubit because the state on final row of 
control qubit entangles with both upstream and down- 
stream part of wavefunction on target qubit. Numerical 
calculation on \tPq) agrees exactly with this analysis. 

The first excited state should be such a state that on 
the first row of one of the two interacting qubits, lipi) is 
orthogonal to |'0o)j and its amplitude is nearly zero due 
to large on-site potential from boundary Hamiltonian /i*^. 
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while on the first row of the other qubit is at the 
same state as |^o) with only amplitude modified. If one 
knows the overall amplitude of wavefunction in upstream 
part on the qubit whose first row state is orthonormal to 

— ~ II IV'l) tipsiream II, (8) 
X 

then the scale of energy gap can be found by same rea- 
soning for single qubit: 

A oc e/x^. (9) 

If IV'i) maintains similar weight distribution to |V'o), we 
find that l^/ii) should be orthogonal to |V'o) on the first 
row of the control qubit, in which 1/a; cx while on 

target qubit oc 1/A, thus the energy gap A oc e/A^. 
The larger value l/x on target qubit doesn't determine 
the energy scale because there at the first row has 
appreciable amplitude. 

Numerical calculation agrees with above analysis for 
We also found that the second excited state |^/'2) 
has the same energy gap scaling as \ipi). |'02) on the first 
row of target qubit is orthonormal to IV'o) with nearly 
zero amplitude, while its first row state on control qubit 
is the same as |V'o)- 1/a; on both target qubit and control 
qubit at 1-02) are of the same order 1 /A^, hence the energy 
gap is also proportional to e/A^. The parameter 1/a; on 
target qubit in \'ijj2) is different from those in IV'o) and 
IV'i) because IV'2) at the first row of target qubit has to be 
nearly zero, according to h{CNOT), the suppression of 
amplitude of wavefunction requires larger difference over 
two neighboring rows to supply same amount of kinetic 
energy. In both and |'02)i the value l/x maintains 
the same order as ground state on the qubit whose first 
row state is the same as |V'o)- 

If projection Hamiltonians are applied instead of boost 
Hamiltonians, and at is appreciable, where |^) 

is the state on the row before h{B,X), then A oc e/A^ 
due to the same reason as for single qubit. If there is 
one boost Hamiltonian and one projection Hamiltonian 
on the two interacting qubits, energy gap is still A oc 
e/X'^. It is interesting to note that by projecting the 
target qubit into one state on its last row and boost this 
state, one can manipulate the state on the final row of 
the control qubit, and this can be observed in Eq.©. 

Numerical calculations confirm our analysis. For 
example, the energy gap for two 4-row qubits ended 
with h{B,X), as shown in Fig.QJi) with = 10e(/ — 
cr^), /i^ ^ 10e(/ + 0-3), is A 0.0574, 2.43 x 
10-3, 3.05 X 10-^ 3.12 X 10-^ 3.13 x lO^^e at A = 
1, /To, 10, ^/T000, 100 respectively. 

In the works by Mizel et.al.'Tl'l, when h{B,X) are ap- 
plied to two interacting qubits, their lower bound energy 
gap agrees with result in the present paper, however, they 
didn't elaborate further for multiple interacti ng qubits. 
Fig. (7) in |0| shows how to use teleportation(16j to in- 
crease A for single qubit evolution, and there A is pro- 
portional to e/A^, which is incorrect because even two 



interacting qubits will give energy gap of e/X^ as shown 
in our analysis and numerical calculation for two qubit 
case, while the teleportation circuit involves more than 
two qubits interacting with each other. However, contin- 
uous single qubit evolution can always be combined into 
one unitary transformation, there is no need to evaluate 
this situation. 

With multiple qubits interacting with each other, we 
need to evaluate on the top part of each qubit the param- 
eter 1 /x defined in Eq. (jSJ) assuming that only on the first 
row of that qubit the excited state is orthonormal to Itpo) 
while states on all other qubits remains the same as cor- 
responding ground state with only magnitude changed. 
Then the minimum 1/a; gives the energy gap scale as 
e(l/a;)^j„. Thus the energy gap depends on the detail of 
the circuit. 

In principle, when all qubits end with either projection 
or boost Hamiltonian containing same A, at ground state 
and first excited state, the boost Hamiltonian or projec- 
tion Hamiltonian raises the amplitude of wavefunction 
on the final row by A compared without boost or projec- 
tion Hamiltonian. When estimating 1/a; of Itpi) for any 
qubit, say qubit A, the boost Hamiltonian, not the pro- 
jection Hamiltonian, on qubit A itself contributes A to 
X, and projection or boost Hamiltonians from its control 
operation mates contribute A to its x value; contribution 
from those qubits not directly interacting with qubit A 
can be analyzed according to ground state wavefunction 
distribution, Eq.Q. 

For example, in Fig. ^p) qubit 1 CNOT controls qubit 
2, and in its downstream part, qubit 2 CNOT controls 
qubit 3, so on to qubit 5. When all qubits end with 
h{B^ A), we find that, for A 3> 5 ( 5 is the typical length 
of qubit) the parameter in Eq.© on qubit 5 has contri- 
bution only from h{B,Xys on itself and qubit 4, those 
on qubit 3, 2 and 1 doesn't make contribution because 
on qubit 5 its upstream part can entangle with the fi- 
nal rows of other qubits except for qubit 4. We get 
X oc A^, A^, A'', A^A5 for qubit 5, 4, 3, 2 and 1 re- 
spectively. The energy gap is determined by value of 1 /x 
on qubit 1 or 2, hence A oc e/A^". As number of qubits 
N in such a chain increases, the energy gap shrinks as 
e/A2^. 

The energy gap of a GSQC circuit may be exponen- 
tially small depending on detail of circuits. The chain 
in Fig. 1^3) is the most dangerous circuit giving the least 
energy gap, the other circuit that also suffers exponen- 
tial shrink of energy gap A oc e/A^^ is a qubit interacting 
with N other qubit as the control qubit. These two kinds 
of circuits are quite common in quantum algorithm, such 
as quantum Fourier transform. 

In order to break such exponential shrink of A like in 
Fig.(^), the circuit is modified by inserting teleporta- 
tion boxes for all qubits between two control Hamilto- 
nians, as in Fig.©, which shows how the circuit in the 
dotted box of Fig. (^) is modified. A teleportation box 
teleports quantum state of a qubit to the other one, and 
it can be implemented in both standard paradigm|l6j 
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FIG. 2: Modified circuit of the dotted box in Fig.l^), in 
which four teleportation boxes are inserted before and after 
HiCNOT). The circuit inside each teleportation box is sim- 
ilar to the teleportation circuit in [T^ . Label / represents for 
identity transformation Hamiltonian h{I), H for Hadamard 
transformation Hamiltonian h{H) and P(0) for projection 
Hamiltonian h{\Q),\). The dotted lines represent for part 
of circuit not showing. 



and GSQClIi. The whole chain in Fig. lit) should be 
modified similarly. In Fig.©, the four boxes contain tele- 
portation circuits similar to that in 15]. On the middle 
qubit of teleportation boxes in Fig.(j2l, the row labeled 
by M opens a gate for its left side qubit 's upstream part 
to entangle with downstream rows at its right side qubit, 
thus the exponential shrink of 1 /x along the chain is pre- 
vented. The minimum 1/x happens on qubit 3 or qubit 
6 in Fig.|j21l. For example, on qubit 3, the parameter 
1/x has contributions only from projection Hamiltoni- 
ans /i(|0),A) on qubit 1, 2, 6 and 7, hence 1/x cx 
While projection or boost Hamiltonian on qubit 8 or 10 
doesn't contribute to 1/x on qubit 3 because through 
rows marked by AI on the middle qubit in teleporta- 
tion boxes, their final row states entangle with the first 
two- row state of qubit 3, and they surely also entangle 
with the last row state on qubit 3. Qubit 9 doesn't con- 
tribute to 1/x on qubit 6, thus has no effect on qubit 
3 because the boost or projection Hamiltonian on con- 
trol qubit doesn't affect the parameter on target qubit 
at ground state, while the lowest excited state remains 
similar with ground state on qubit 9 if only on qubit 3 
the first row state is orthonormal to |'0o)- Just like qubit 
8 and 10, other qubits further away that do not show up 



in Fig. 121) do not affect 1/x on qubit 3. By inserting tele- 
portation boxes, we find that when the amplifying factor 
of all boost and projection Hamiltonians is A, the energy 
gap for arbitrary GSQC circuit is always 

A cx e/A^ (10) 

Hence there is no chain action on 1/x like in Fig(^) 
that leads to exponential shrink of energy gap. In fact, 
the upper right and the lower left teleportation boxes are 
not needed concerning on circuit like Fig.(^3), however, 
I keep them for more general consideration. 

For arbitrary circuit, it can always be modified similar 
to Fig.(|2I) and the energy gap will be kept at the scale 
of e/A^. The price of this modification is that the to- 
tal number of qubits increases, thus from measurement 
concern, it requires larger value of A to make sure appre- 
ciable probability to find all electrons on the final row's of 
all qubits. However, this price is worthy because it may 
change energy gap from exponentially small to polyno- 
mially small, while the value of A is only polynomially 
increased. This can be demonstrated in the following 
example. 

Now we can check the energy gap when a quantum 
algorithm is implemented by GSQC. The most powerful 
quantum algorithm up to now is quantum Fourier trans- 
form, which is used in factorization, period finding, etc. 
The detail of quantum Fourier transform, more precisely, 
the inverse quantum Fourier transform, can be found in 
many literatures, such as |0. It takes 0(7V^) control 
operations h(CR\.) to carry out inverse quantum Fourier 
transform by standard time dependent approach, where 
N is the number of qubits involved in the problem, and 
is the unitary transformation shifting phase of state 
|l)by-27ri/2'='tl6i]. h{CRl) can be obtained by replacing 
NOT operator by R^^. in Eq. . One can simply replace 
the time evolution by circuit array and calculate the en- 
ergy gap of the circuit. By inserting teleportation boxes 
between control operations, when all qubits end either by 
/i(|0), A) or h{B, A), the energy gap is A oc e/A^. Because 
after experiencing a control operation each qubit is tele- 
ported to a new qubit, the total number of qubits is of 
the order iV^ instead of N . With each qubit ended with 
boost or projection Hamiltonian, at the probability 
of finding electron on the final row on any qubit is larger 
than (1 — C/A^) with C < 8 being the number of rows on 
the qubit. In order to find all electrons on final rows with 
appreciable probability P > (1 - C/\^Y^ , it requires 
A ~ ^/DN, here we assume there are totally FN'^ qubits 
involved, and D is a constant that could be tuned. Thus 
the probability is P w e^^'-^/^ as N being large number, 
the energy gap is A c>c e/ (D'^N^), noting that standard 
paradigm needs only 0{N'^) steps. 

We conclude that the energy gap of ground state quan- 
tum computer is determined by the number of control 
operation, and any quantum algorithm implemented by 
standard paradigm can be implemented by GSQC, whose 
energy gap is A cx l/iV""^, where N is the number of bits 
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in the problem, and N'^ is the number of control opera- 
tions needed in standard quantum paradigm. There are 
various approaches to make GSQC stay at it's ground 
state, and the time cost is determined by different ap- 
proaches. Estimation of time can be obtained by adi- 
abatic approach, in which we increase A slowly from 1 
to 0(A^*^/^) in all boost and projection Hamiltonians on 
final row of all qubits, GSQC will stay at ground state 



and time cost is of the order N^'' , determined by inverse 
of square of minimum energy gap. By local adiabatic 
approach, the time cost is even shorter^. 
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